Coupling constant constraints in a nonminimally coupled phantom cosmology 
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In the paper we investigate observational constraints on coupling to gravity constant parame- 
ter £ using distant supernovae SNIa data, baryon oscillation peak (BOP), the cosmic microwave 
background radiation (CMBR) shift parameter, and H(z) data set. We estimate the value of this 
parameter to constrain the extended quintessence models with nonminimally coupled to gravity 
phantom scalar field. The combined analysis of observational data favors a value of £ which lies 
in close neighborhood of the conformal coupling. While our estimations are model dependent they 
give rise to a indirect bound on the Equivalence Principle. 
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At the present, scalar fields play the crucial role in 
modern cosmology. In the inflationary scenario they gen- 
erate an exponential rate of evolution of the universe as 
well as a density fluctuations due to vacuum energy. 

Observations of distant supernovae support the cos- 
mological constant term. But two problems emerge in 
this context. Namely, the fine tuning and cosmic coin- 
cidence problems. The lack of some fundamental mech- 
anism which sets the cosmological constant almost zero 
is called the cosmological constant problem. The second 
problem called "cosmic conundrum" is the question why 
the energy densities of both dark energy and dark mat- 
ter are nearly equal at the present epoch. One of the 
solutions to this problem offers the idea of quintessence 
0, 0] which is a version of the time varying cosmological 
constant conception. 

All these models base on the assumption that there is 
a minimal coupling of scalar field to gravity (£ = 0). 
This a priori assumption requires some justification. 
There are many theoretical arguments suggesting that 
the non-minimal coupling (NC for short) should be con- 
sidered. The nonzero £ comes from quantum corrections 
Q, renormalization of classical theory that shifts it to 
one with nonzero £ [1] , in relativity the value of £ = 1/6 
(conformal coupling) is distinguished. Only in this case 
the Einstein equivalence principle is not violated (for de- 
tails see [!])• 

It is expected also that the value of coupling constant 
£ should be fixed by the physics of the problem and there 
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is no free parameters, and then we obtain answer to the 
question "what is the value of £?" Q. However, the 
answer to this question differs according to the theory 
of the scalar field employed. For example, in Einstein 
gravity with the polynomial quartic potential function 
and back-reaction, the value of NC was found [?J or a 
Higgs field in the standard model gives value of NC non- 
positive or larger (or equal) than 1/6 [8(. 

The main goal of this paper is a statistical estimation 
of the coupling parameter from the astronomical obser- 
vations. For this aim we consider the spatially flat FRW 
model where the source of gravity is a noninteracting 
mixture of dust matter and non-minimally coupled phan- 
tom scalar field. 

We treat £ as a free parameter in the model and we are 
looking for constraints on it's value from observational 
cosmology. For simplicity it is assumed a simple and 
natural form of the potential function of the scalar field. 
For the vanishing coupling constant (minimal coupling) 
this potential corresponds to chaotic inflation [9j. This 
paradigm of inflation can be extended by including NC 
[6j . Also the scalar field non-minimally coupled to gravity 
are simple, non-exotic models of phantoms give rise to 
superacceleration (loj . 

Tsujikawa and Gumjudpai TlJ also investigated con- 
straints on the NC parameter from the CMB observa- 
tions. Recently Jankiewicz and Kephart fl2| applied the 
method of WKB approximation to study the NC in the 
FRW universe with the scalar field and found the best 
fit value the NC parameter very close to 1/6. Futamase 
and Maeda [l3[ in the case of the quadratic potential of 
the inflaton field found that modulus of NC less than 
10~ 3 leads to sufficient inflation (see also [l4[ for other 
constraints). 

In order to determine the NC parameter we use a re- 
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cent expansion history of the universe rather than early 
stages of its evolution (inflation). In estimation of the pa- 
rameter £ we do not assume any arbitrary parametriza- 
tion for the equation of state parameter w(z), but we 
derive its exact form directly from the model (dynam- 
ics). While this equation has in general very complex 
form, in practice its simple approximations are used. If 
we expand the w{z) relation in the Taylor series with the 
respect to the redshift z then we recover the well known 
parametrization in the linear form of the redshift or the 
scale factor. All the coefficients in the series depend on 
the parameter £ and initial conditions. Then after con- 
straining initial conditions from the observational data 
we can obtain corresponding limits on £. 

We assume the flat model with the FRW geometry 
with the signature (—,+,+,+) and a source of gravity 
is the phantom scalar field ip with an arbitrary coupling 
constant £. The dynamics is governed by the action 
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where m 2 p 



(1) 

(8irG) 1 ; for simplicity and without lost of 
generality we assume AttG/3 = 1 and m 2 = 1. 

In our paper [l5[ we studied generic features of the 
evolutionary paths of the flat FRW model with phan- 
tom scalar field non-minimally coupled to gravity. We 
reduced dynamics of the model to the simple case of au- 
tonomous dynamical system on the invariant submani- 
fold (tp, ij) f ) (a prime denotes differentiation with respect 
to the natural logarithm of the scale factor) and we 
founded that principally there is one asymptotic state, 
which corresponds to the critical point in the phase space 
i/jq = il/v 7 ^ and ip' =Q. This critical point is also the 
de Sitter state (w = —1). There are two types of scenario 
leading to this Lambda state (depending on the value of 
£), through 

1. the monotonic evolution toward the critical point 
of a node type for < £ < 3/25; 

2. the damping oscillations around the critical point 
of a focus type for 3/25 < £ < 1/3. 

In both evolutionary scenarios we obtained linearized 
solutions of the dynamical system in the vicinity of the 
critical point. Using relation of the equation of state pa- 
rameter on the dynamical variables ip and ip 1 we were able 
to derive exact forms of w(z) relations in both scenarios 
(see [HI for details). 

Now, we present method of estimation model param- 
eters from combined analysis of the parameters wq, u>i 
and W2 from the combined astronomical data for three 
parametrization of w(z) (see Table HJ. 

To constrain the unknown values of model parame- 
ters we used the the set of N\ = 192 SNIa data, CMBR 
shift parameter, measurement of the parameter A coming 
from the SDSS luminous red galaxies and N 2 = 9 obser- 
vational H(z) data. The likelihood function for such data 



sets has the following form 



C{9) tx exp 



(2) 



where 9 = (Ho,0, m o,w ,Wi,w 2 ) and % 2 (6>) = Xsn(®) + 

X 2 r (O) + X 2 a(0) + X 2 h(O)- 

Let us consider quantities from the definition of func- 
tion x 2 - Using the SNIa data [H, [13, [3 we base on 
the standard relation between the apparent magnitude 
m and luminosity distance d^: m — M = 51og 10 D^+Ai, 
where M is the absolute magnitude of SNIa, Ai = 
— 51og 10 i7 + 25 and Dl = HqcIl. The luminosity 
distance depends on the cosmological model considered 
and with the assumption that k — is given by g?l = 
(1 + z)c f Q h(z>) ■ ^he X%h function has the following 
form 



xh(0) 



E 

i=l 



obs 



th 



(3) 



where n° bs = m, - M, [if = 5 log 10 D Ll + M. 

As we have written before we also used constraints 
which comes from the so called CMBR shift parameter 
which is defined as R th = y^K^ f* doc jffadz. The x\ 
has the following form 



X 2 R (0) = 



R obs - R 



Ih 



(4) 



where R obs = 1.70 ± 0.03 for z dcc = 1089 (HH^. 

We also used constrains coming from the SDSS lumi- 
nous red galaxies measurement of A parameter (A obs = 
0.469 ± 0.017 for z A = 0.35) [U, which is related to the 
baryon acoustic oscillation peak and defined in the fol- 

lowing way A th = ^0 (^) "* [£ C ^dz ' 
The x\ function has the following form 



X\(0) = 



A 



iii 



A 



obs 



(5) 



Finally we add constraints coming from observational 
H(z) data (N = 9) HHHIH. This data based on the 
differential ages (dt/dz) of the passively evolving galaxies 
which allow to estimate the relation H(z) = a/a = — (1 + 
z)~ 1 dz/dt. The function Xh nas ^ ne following form 



X 2 H (0) = 



N 2 

E 

i=l 



H(zi) - Hi(zi) 



(6) 



After the marginalization likelihood function over the 
parameter Hq in the range < 60, 80 > the values of 
model parameters (Ct mt o,wo,Wi,Wi) were obtained via 
the x 2 minimization procedure and are gathered in Table 
UJ There are also the la uncertainties which correspond 
to the parameter value where the likelihood function dif- 
fer from the best fit likelihood by a factor of e = . 
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TABLE I: Estimations of expansion coefficients for three different parametrization of w(z). 



Case 


parametrization 


Wo 


Wl 


W2 


(1) 


w(z) — Wo ± W\Z ± W2Z 


-1.29 ±0.18 


1.50 ±1.03 


-0.57 ±0.75 


(2) 


w(z) = w + w 1T ^ + w 2 (j^) 2 


-1.38 ±0.13 


2.75 ±0.87 


-2.13 ± 1.33 


(3) 


w(z) = wo + Wi In (f + z) ± W2QXL (1 ± z)) 2 


-1.34 ±0.24 


2.16 ± 1.71 


-1.25 ± 1.73 



TABLE II: Solutions to the system of equations w x , w x and 
w\ with values of wo, mi and W2 from Table [I] 



Case 


Xo 


yo 


yo/xo 




(1) 


7.45 ± 7.05 


-8.55 ± 7.77 


-1.15 ± 2.13 


0.14 ±0.03 


(2) 


6.44 ± 2.75 


-8.22 ± 3.11 


-1.28 ± 1.03 


0.15 ±0.02 


(3) 


6.53 ±8.77 


-8.09 ± 10.12 


-1.24 ±3.21 


0.15 ±0.05 



For every fitting parametrization we obtain a set of 
three numbers (wo, w\,W2) (see Table UJ). Then basing on 
expansion in the Taylor series of the exact forms of EoS 
functions up to second order term (w x , w x , w x ), we can 
solve the system of three equation for three values £, xq 
i yo- The errors of Axq, Ayo and A£ are calculated in a 
standard way. From Table [TT] we can see that the param- 
eter £ is closed to the its conformal value 1/6. Therefore 
the combined analysis of observational data favors a value 
of £ which lies in close neighborhood of the conformal 
coupling. The second conclusion is that the oscillatory 
scenario is favored. In the case of parametrization (2) we 
can exclude the monotonic scenario on the la confidence 
level. 

Finally we consider parameters wq, iui, W2 as functions 
of parameters £, Xq, yo- The values of model parameters 
^m,0j xq, yo obtained via the \ 2 minimization proce- 
dure together with la errors are gathered in Table IILTI 

The initial condition (&o,yo) defines the initial rate of 
change of scalar field (ln(V>j — V^o))' = ifr'Jtyi — ipo) = 
Vq/ x o = (Inxo)'- This value is negative and is presented 
in the third column of Table [TT] and in the fourth column 
of Table [TTTJ 

We can see that in this case the values of parameters do 
not differ from the values obtained by indirect method via 
the estimation of wq, w\, W2 but here the uncertainties of 
the parameter values are smaller. In Table HVl we gath- 
ered the values of the parameter of non-minimal coupling 
£ estimated via direct method together with lcr, 2a and 
3a confidence levels. The value of £ = 1/6 is consistent 
with our results at a 3a confidence level for parametriza- 
tion (1) and at a 2a confidence level for parametrization 
(2) and (3). As one can conclude the oscillatory scenario 
is favored, moreover, we can exclude the monotonic sce- 
nario at a 2a confidence level for the parametrization (1) 
and (3) and at a 3a confidence level for parametrization 
(2). 

While the standard quintessence idea is most popu- 
lar one, scalar field cosmology with minimally coupled 



scalar field require existence of a special conditions for a 
tracker solution and potential functions of the scalar field 
[25l | . The standard quintessence evolutionary scenario 
has been extended by introduction of the non-minimal 
coupling between gravitation and scalar field. Hence the 
cosmology with such a fields contains an additional pa- 
rameter which should be determined from a observational 
data. We demonstrate that astronomical data and cos- 
mography may be very useful in estimation of this free 
parameter of a model. 

In our approach the effects of non vanishing coupling 
constant are dynamically equivalent to the effects of sub- 
stantial dark energy characterized by the coefficient of 
the equation of state. The dynamics with positive value 
of the NC parameter admits two types of evolutionary 
scenarios: the monotonic evolution and damping oscilla- 
tions around the stationary de Sitter state. Given the 
estimated value of the parameter £ we found that the 
oscillating scenario is favored. 

We expanded the exact form of the EoS parameter up 
to second order term and then constrained the effective 
form of the EoS which determines directly the value of 
the NC parameter. As a result performed combined anal- 
ysis we obtain that value of £ from neighborhood of con- 
formal coupling is favored by observational astronomical 
data, i.e. SN la data, Baryon Oscillation Peak, CMBR 
shift and H(z) data. It was demonstrated that value of 
estimated NC parameter does not depend on used expan- 
sion formula up to the second order although errors can 
be different. 

The General Relativity in based on the Einstein Equiv- 
alence Principle, which includes Weak Equivalence Prin- 
ciple, Local Lorentz Invariance and Local Position In- 
variance. Therefore if this universality condition is valid 
then gravity can be described in terms of Ricmannian ge- 
ometry. Our estimation of coupling constant is based on 
the FRW model filled by non-minimally coupled phantom 
scalar field but one should also note the existence of the 
direct constraints on WEP. Then deviation value of cou- 
pling constant from the conformal coupling can measure 
the deviation from the WEP and the GRG in the Einstein 
formulation. Our estimation, based on a model with the 
non-minimally coupled phantom scalar field, gives the 
deviation from the WEP A£/£ ~ 1CT 1 . On the other- 
hand the string theory suggest the existence of strength 
scalar field (dilations) whose couplings to matter can vi- 
olate the WEP 26]. This provides a new motivation for 
high-precision experiments verifying the universality of 
free fall to the KT 12 level [H [21, [29[. The lunar Laser 



4 



TABLE III: The values of directly estimated parameters Q m ,o, f, xo, yo after substituting w x , w x , w x into the H(z) formula 
for three different parametrization of w(z). 



Case 


"m,0 


Xq 


Vo 


yo/xo 




(1) 


0.27 ±0.02 


7.33 ± 0.90 


-8.43 ±0.75 


-1.15 ±0.24 


0.14 ±0.01 


(2) 


0.27 ±0.02 


6.41 ± 0.32 


-8.18 ±0.32 


-1.28 ±0.11 


0.15 ±0.01 


(3) 


0.27 ±0.02 


6.60 ± 0.72 


-8.16 ±0.81 


-1.24 ±0.26 


0.15 ±0.01 



TABLE IV: The values of the parameter of non-minimal cou- 
pling £ estimated via direct method together with la, 2a and 
3a confidence levels. In the cases (1) and (3) we can exclude 
the monotonic scenario at a 2a and in the case (2) at a 3a 
level. 



Case 




la 


2a 


3a 


(1) 


0.14 


0.01 


0.02 


0.04 


(2) 


0.15 


0.01 


0.02 


0.03 


(3) 


0.15 


0.01 


0.02 


0.05 



Ranging experiment has also verified that the Moon and 
the Earth fall with the same acceleration toward the Sun 
to better than one part in 10 12 precision [30]. Also ob- 
served limits on evolution of the binary pulsar B0655±64 
orbit provide new bounds on the violation of SEP [3lj | 
(for the description of current and future projects for the 
improvement of the accuracy of the experiments as com- 
pared experiments on ground see [3^]). 

This work has been supported by the Marie Curie Ac- 
tions Transfer of Knowledge project COCOS (contract 
MTKD-CT-2004-517186). 
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